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Abstract — In this work, the boundary value problem of 
simply supported rectangular Kirchhoff plates subjected to 
applied transverse loads is solved by the method of finite 
Fourier sine transform. The finite Fourier sine transform 
method was adopted as the analytical research tool due to 
the Dirchlet boundary conditions of the plate problem. 
Application of the finite Fourier sine transform to the fourth 
order governing partial differential equation of the 
Kirchhoff plate problem and the associated boundary 
conditions simplified the problem to an algebraic problem in 
the transform domain. The solution is obtained in the plate 
domain by inversion. The problem was solved for general 
distributed load p(x, y), point load applied at an arbitrary 
point on the plate, uniformly distributed patch load over the 
plate region xo< x <xi, yoNy <yi, and uniformly distributed 
load over the entire plate. The finite Fourier sine transform 
solutions obtained in each case were found to be identical 
solutions obtained with the Navier’s double trigonometrical 
series method as presented in Timoshenko and Woinowsky- 
Krieger. The finite Fourier sine transform method was found 
to yield exact solutions to the classical thin plate flexure 
problem for simply supported edges. 

Keywords — Finite Fourier sine transform method, 
Kirchhoff plate, Dirichlet boundary conditions, distributed 
transverse load, patch load, point load, Navier’s double 
trigonometric series method. 

I. INTRODUCTION 

Plates are three dimensional structural members with 
extensive applications in civil, mechanical, aeronautical, 
naval and geotechnical engineering used to carry external 
loads by the development of bending resistance about the 
two axes of the plate [1, 2, 3], 

The term plate theory denotes an approximate theory used to 
determine the stress fields and deformation field in elastic 
bodies one dimension of which (the plate thickness, h) is 
small compared with the other dimensions (the width and 
length of a rectangular plate supported at the edges or the 
diameter of a circular plate) [4] . The approximations consist 
of the introduction of certain simplifying assumptions into 
the governing kinematic, stress strain and equilibrium 


equations of the mathematical theory of elasticity [3, 4, 5], 
These simplifications yield results which do not differ 
significantly from those obtained from the exact equations 
for the range of definition of the problem. The 
simplifications used in various plate theories derive from the 
definition of a plate as a three dimensional structure with 
one small dimension; and also from the consequences of 
Bernoulli-Navier’s hypothesis for beams when extended to 
plates. 

In the classical Kirchhoff -Love’s plate theory, the influence 
of transverse shear strains is assumed to be negligible, and a 
simultaneous consideration of kinematics, stress-strain law 
and the differential equation of equilibrium for an 
infinitesimal plate element results in a fourth order partial 
differential equation as the governing equation of 
equilibrium [7], Consequently, the number of boundary 
conditions appurtenant to the support conditions appears to 
be in disagreement with the order of the governing partial 
differential equation [3, 8], This limits the validity of the 
expressions for the shearing forces to the open region of the 
plate middle surface and introduces Kirchhoff s shearing 
forces for the boundary of the plate. Three actual boundary 
conditions at each edge of the plate have to be replaced by 
two approximate conditions transformed in the Kirchhoff 
sense [2, 4], 

Despite the shortcomings of the classical Kirchhoff-Love 
plate theory, it is well documented that for the majority of 
engineering applications, the theory gives sufficiently 
accurate results. The limitations and imperfections of the 
classical Kirchhoff-Love plate theory have led to the 
development of other plate theories. Some of these are 
Reissner plate theory [9, 10]; Mindlin plate theory [11], 
Henky refined plate theory [12, 13], Shimpi refined plate 
theory [14], Higher Order Plate Deformation theory [2], 
Third Order plate theory [15], Leung’s Plate theory and 
Osadebe plate model [16]. Modified plate theories have also 
been used in plate bending analysis [17]. 

The plate problem in general is a boundary value problem 
which is a system of differential equations to be satisfied in 
the plate domain and the associated boundary conditions to 
be satisfied at the plate boundaries [18], The plate problem 
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has been solved successfully in the technical literature using 
two basic approaches - classical methods - Navier’s double 
trigonometric series and Levy’s single trigonometric series 
methods, and Numerical or Approximate Methods - Finite 
Difference Method [17], Finite Element Method, Boundary 
Integral Method, Variational Methods [20, 21] (Ritz 
Variational method [22], Galerkin Variational Method), 
Integral transform methods (Laplace transforms, Fourier 
transforms and Hankel transforms) [23] and Perturbation 
methods. In this work, the finite Fourier sine transform 
method is applied to solve the boundary value problem of 
simply supported rectangular Kirchhoff plates under given 
transverse loads. 


II. METHODOLOGY 

The finite transforms follow from the theory of Fourier 
series [24, 25] 

The finite sine transform S„ of a function of xf(x) is defined 
as 

n l 

n _ . „ , . . Z f . . , . TITZJC - 

S n = S(f(x)) = - f(x)sm—dx ; ( 1 ) 

l J I 

o 

m= 1, 2, 3 ..., n = 1, 2, 3 ... 

where 0 < x < /, S n = S(f(x)) is the finite Fourier sine 

transform of fix) 

with its inverse sine transform as 


S-\S n ) = f(x) = Y j S n sin^ (2) 

17=1 


where S 1 is the inverse finite sine transform 
The finite sine transform is commonly used with Dirichlet 
boundary conditions, that specify the value of f(x) at the 
domain boundaries;,!: = 0, and x = l. 

The Finite cosine transform C„ of a function of xf(x), 
normally used with Neuman boundary conditions that 


df 

specify the value of — at the domain boundaries,! = 0, and 

dx 


x = l is defined as 


9 f 1 IITLX 

C n = C(/(x)) = ~^f(x)cos—dx , (3) 

n = 1, 2, 3 ... 


The inverse cosine transform, denoted by C \ is defined 
by: 


^ oo 

c~\c n ) = f(x) = -j- + COS ^ (4) 

77=1 

These transforms reduce a partial differential equation PDE 
to an ordinary differential equation ODE. 


If fix, y) is a function of two independent variables x and y, 
defined in a given region 0 <x<a , 0 <y<b, its double finite 
Fourier sine transform F s (/( x, y)) is defined by [26, 27] 


a b 

F s (f(x, ?)) = F(m, n ) = J J fix, y)sin 
0 0 


mnx . nny , , 

sin — — dxdy 

a b 


...(5) 

The inverse double finite Fourier sine transform is given by 
the double series: 


0° °° / x 

Ff'(f(m,n)) = f(x,y) = — ^^F(m,n)sin^^^Jsin^, 

m n 

m, n = 1, 2, 3 •••(6) 

where F s ' ( f (m, n)) is the inverse double finite Fourier 

sine transform of (f(m, n )) 


III. APPLICATION OL THE LINITE FOURIER 
SINE TRANSFORM METHOD 

The governing partial differential equation to be solved is 
given by Kirchhoff plate equation: 

3 4 w d A w d 4 w p(x,y) 

— r- + 2— — - + — - ^ ‘ v (7) 

Sr 4 dx 2 dy dy D 

where p(x, y) is the transverse distributed load acting on the 
plate, D is the flexural rigidity of the plate material. 

Eh 3 

D — — ,E = Young’s modulus of elasticity, h = 

12(1 -ft 2 ) 

plate thickness, p is the Poisson’s ratio, and vv(x, y) is the 

transverse deflection of plate, x and y are the space variables 

in the plane of the plate, andO <x<a, 0<y<b. 

The boundary conditions for simple supports at the plate 

edges: x = 0, a; y = 0, b are w = 0on.v= 0, a,y = 0, b, w xv = 

0 on x = 0, a , w yy = 0 on y = 0, b 

d 2 w d 2 w 

where w xx = — y-, and w yv = — - 

ox" ' ' dy 


Transverse Distributed load p(x,y ) 

For distributed transverse load p(x, y), taking the finite 
Fourier Sine transforms of both sides of Equation (7), we 
have Equation (8) 


b a 

II 

0 0 


0 4 w 

dx 4 


+ 2 


d 4 w 
dx 2 dy 2 


+ ■ 


d 4 w ^ 


dy 


sin a m x sin P n y dxdy 


D 


II 


p(x, y) sin a m x sin P n y dxdy 


( 8 ) 


where Ct m 


mn 

a 


P n 


nn 

b 
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Using the linearity property of the finite Fourier Sine 
transform, and noting the finite Fourier Sine transform of 

— — given by Equation (9) 
dx 4 


0 0 


1 1 — sin a m x sin P n y dxdy 
b 

| [(- 1 )" 1 w xr 1 x=a ~™xx 1 *=0 ] sin fitly dy 
0 

a b 

J J w(x, y) sin a m x sin (i n y dxdy 
• ••(9) 


+a„ 


o o 


= a m w(tn, n ) 

since vi ’ , , I = WvaL-o = 0 

where 

b a 


w(m, n) - | | w(x, y)sin a,„xsin P n y dxdy (10) 

a 0 

w(m, n ) = finite Fourier Sine transform of the transverse 
deflection function w(x, y), we have Equation (11): 


b a 

1 1 (a 4 + 2 af„P^ +p 4 )w(x, y)sin a m xsin P H y dxdy 
0 0 

1 bo 

= — p(x, y) sin a m x sin P n y dxdy ( 11 ) 

U j j 

a 0 

2 b a 

( a /n + Pn) J J w(x, y) sin a,„x sin P„y dxdy 
o o 

1 b a 

- — 1| p(x, y)sina m xsinP, ; y dxdy ( 12 ) 


0 0 


( a m + Pn) 


D 


(13) 


where p m „ is the finite Fourier Sine transform of the 

transverse distributed load 

Hence, 

Pmn 


W 


mn D( a 2 m +V 2 n ) 2 
By inversion. 


(14) 


w(x, >0 - — — V V w mn sin a m xsin p„y (15) 
Dab , , 

m = 1 n = 1 

w(x, y ) = — y y p,nn sma, » xsi ” fi n y (16) 

' abD ^n 


. mux . nny 
4 ® “ Pmn sin sm - 

zz a b 


abD ^ 4 ~ l f 

m n 


f mn\ 
VV a J 


,2\ 2 

nn | 


(17) 


n’abD 


. mux . nny 

oo oo p sin sin — — 

y^ y^ a b_ 

L—t / 2 2 ^ 

m n 

— 2 + 7 T 

\a b J 


(18) 


This satisfies the boundary conditions 

Case of concentrated load p(x, y) = Pq at x = q, y = r\ 
inside the plate domain 

The load p(x, y ) is represented using Dirac delta functions as 
p(x, y) = P 0 5(x - £,)S(v - r|) (19) 

where Po is a constant and S(x — cj8(y — r|) are the Dirac 
delta functions, \ and r) are the coordinates of application 
of point load on the plate, 
then, 0 < (; < a; 0 < r| < b 
b a 

1 1 V 4 w(x, y) sin a m x sin P n y dxdy 

o o 

b a 

= — 1| Po&(x ~ £)8(y - Tl) sin a,„xsin P„y dxdy (20) 


0 0 


((X„, + 2(X m P n + fin ) w n 


Pq . mnt . mrn 

— sm sm 

Dab 


Pq . mn't . 77701 

— sin sin 1 

Dab 


s \2 

\ 2 

f 77771 I 

nn i 

+ 


A a J 

,bj J 


(21) 


( 22 ) 


Dn 


. Hint. . 77701 

sin sin 


4 f 2 l\ 

777 77 


(23) 


By inversion, 

^ GO GO 

W(X, y) = — ZZ W '„ n Sin a r,r XSm P,A 
m n 

m,n = 1 , 2 , 3, 4 . . . 


(24) 


. 7777I& . 77701 . 7777IX . 777tV 

. oo co „ sm sm 1 sm sm — — 

4 ^ ^ ‘o a b a b 


ab Dn 

m n 


^2 2\ 

777 77 

2 + 77 
\a b 


...(25) 
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4Po 

Dabn 4 


00 00 sin 

zz- 

m n 


mnt . hkt\ . mux . tiny 

sin 1 sin sin — — 

a b a b 


/ 2 2 
m n 

— 2 + 72 
V a b J 


\2 


...(26) 

w(x, y ) = P 0 K(x, y, t, n ) (27) 

K = Kernel or Green function. This, satisfies the boundary 
conditions. 

The bending moments are found from the bending moment 
displacement relations. 

For point load P at the centre of the plate, 
x — t, = a! 2, y = r\ — b/2 the deflection, and loading 


moment value become 



m= 1,3,5 . . . 


4^o 

Dn^ab 


. 9 mn . 9 mi 
00 00 sin' sin — 


ZZ 

m n 



2 


b 2 J 
...(28) 


n = 1, 3, 5 ... 


M v 


M 


yv 


4Pa 

~1F 

4Pa 

IN 


ZZ 

m n 

00 00 

zz 


t 1 \tn~vn — 2/ 2 2 . r, 2\ 

(-1) (m a +0.3 n ) 

r 2 2 , 272 

(m a +n ) 

/ 1 — 2/ 2 1/99 2 2\ 

(-1) (n + 0.3m a ) 

r 2 2 , 272 

(m a + n ) 


(29) 

(30) 


Uniformly distributed patch load over the plate region 

xn<jc<x 1 yo<y<yi 

The plate deflection due to a uniformly distributed patch 
load over the region xo<x<xi, yo<y<y 1 is found by integration 
of the point load solution; as 


w(x, y) = 
4 


, . mnt . nnr\ . mux . mny „ , 

00 00 yi x i /?(£,, t]) sin z sin L sin sin dc^dr \ 

' a b a b 


Dabn q 


no- 


2 2 
m n 


_ . mux . nny 

16,, °ZZ " S " TS "4 


n 6 D 


mn 


f 2 

m 


\ a 


2 7 
n 

77 

b y 


_ . mnt . mnu . mir\ . nnv 

= sin sin sin sin 

IIIIL r r\ 1 

a la b 2b 


...(31) 

(32) 

(33) 


where, u = (x 0 - x, ), v = (y, - y 0 ) 

From the bending moment displacement equations, 


c 

1 £ 00 00 inn 

«„=^zz- 

7T 

' v mi m 


(2 2 7 
m n 

2 + ^72^ 

\a b 


. mnx . nny 

sin sin — — 

a b 


mn 


f 2 2 A 2 

m n 

+ 77 

\a b J 


Myy 


16 Po 


c 

go co ^ mn 

zz- 


( 2 
n 


m 




1 2 + 1 1 2 

V b a J 


...(34) 


. mnx . nny 

sin sin — — 

a b 


mn 


(2 2 7 
m n 

~7T + 72 
V a b J 


...(35) 

Uniformly Distributed load po over the entire plate 

The plate deflection for uniformly distributed load po on the 
entire plate domain is found by integrating the point load 
solution over the entire plate are as: 
w(x, y) = 

, . mnt . nnr\ . mux . mny „ . 

00 00 0 a pq sin sin sin sin dtjdx \ (36) 

' a b a b 


Dabn 


-ZZJJ- 

m n 0 0 


m 2 n 2 
+ ? 


. [ mnx i . [ mry | 

. , oo oo sin sin — - 

l6p ° YjYj 

Dtc ™ M 


mn 


(2 2 \ 

m n 

— 2 + 77 
V a b 


01 ) 


m= 1, 3, 5 ... n = 1, 3, 5 ... 

From the bending moment-displacement equations, the 
bending moment distributions become 


m n 

00 CO I 9 f F 


1 S' W UJ 

«,,=^zz- 


2 7 


Dn? 


m n 


t 2 

b 


. mnx . nny 

sin sin — — 

a b 


mn 


( 2 2 A 2 

m n 

~^2 + 77 
V a b J 


...(38) 


, ^ 00 00 

«»=1^ZZ 


f 2 2 7 

n m 

77 + l l ^2 
\b ay 


yy 


DtC 


m, n = 1, 3, 5, 7 ... 


. mnx . nny 

sin sin — — 

a b 


mn 


(2 2 A 2 

m n 

~2 + 77 
\a b J 


...(39) 


IY. RESULTS AND DISCUSSIONS 

The Finite Sine transform method has been applied to the 
boundary value problem of simply supported rectangular 
Kirchhoff plates under general distributed load p(x, y), point 
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load P at ( C, r| ) uniformly distributed patch load over a 
given area of the plate and uniformly distributed load over 
the entire plate. The deflection functions computed are 
shown in each case as Equations! 18), (26), (32) and (37). 
The deflection functions were found to be identical with the 
Navier’s double trigonometric series solution for the 
deflection in each case. The double series of infinite terms 
obtained for the deflection w(x, y) for point loads is a rapidly 
convergent series and the deflection at any point can be 
obtained with good accuracy by considering only the first 
few terms. For point load applied at the plate centre, the 
deflection was observed to be symmetrical about the plate 
axes of symmetry, and maximum deflection was found to 
occur at the plate centre. The series solutions obtained for 
both w(x, y). M x and M y (for point load) were found to 
diverge at the point of application of the point load. 


[Vol-4, Issue-3, Mar- 2017] 
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Variation of maximum deflection at the plate centre, 
with the plate aspect ratio are tabulated inTable 3 

For the case of uniformly distributed load on the entire plate, 
the displacement was found to be identical to the 
corresponding Navier’s double trigonometric series solution. 
The deflection was found to be a double series of infinite 
terms, symmetrical about the two axes of symmetry of the 
plate. Bending moment functions found from the moment 
displacement relations were similarly found to be doubly 
symmetrical about the plate axes. Maximum values of 
deflection and bending moments for various b/a values were 
found to occur at the plate centre and are tabulated as shown 
in Table 1. The convergence characteristics of the series of 
w, Mxx, Myy are shown in Table 2, which shows the 
deflections converge faster than the bending moments. 


Table. 1: Deflection and Bending Moment coefficients for Simply Supported Rectangular Kirchhoff Plates under uniformly 

distributed loads 


b/a 

4 

17 P a 

= F 

Timoshenko and 

W oino wsky-Krieger 

Present study 

(deflection 

coefficient) 

M X x max 

Timoshenko 

and 

Woinowsky- 

Krieger 

Mxx max 

Present 

study 

(moment 

coefficient) 

Myy max 

Timoshenko 

and 

Woino wsky- 
Krieger 

Myy max 

Present 

study 

1.0 

4.06 x 10“ 3 

4.062 x 10“ 3 

0.0479 

0.047886 

0.0479 

0.047886 

1.1 

4.85 x 10“ 3 

4.85 x 10“ 3 

0.0554 

0.0554 

0.0493 

0.0493 

1.2 

5.64 x 10“ 3 

5.64 x 10“ 3 

0.0627 

0.0627 

0.0501 

0.0501 

1.3 

6.83 x 10“ 3 

6.83 x 10“ 3 

0.0694 

0.0694 

0.0503 

0.0503 

1.4 

7.05 x 10“ 3 

7.05 x 10“ 3 

0.0755 

0.0755 

0.0502 

0.0502 

1.5 

7.72 x 10“ 3 

7.724x 10‘ 3 

0.0812 

0.08116 

0.0498 

0.049843 

1.6 

8.30 x 10“ 3 

8.30 x 10“ 3 

0.0862 

0.0862 

0.0492 

0.0492 

1.7 

8.83 x 10“ 3 

8.83 x 10“ 3 

0.0908 

0.0908 

0.0486 

0.0486 

1.8 

9.31 x 10“ 3 

9.31 x 10“ 3 

0.0948 

0.0948 

0.0479 

0.0479 

1.9 

9.74 x 10“ 3 

9.74 x 10“ 3 

0.0985 

0.0985 

0.0471 

0.0471 

2 

10.13 x 10“ 3 

10.12866 x 10 -3 

0.1017 

0.101683 

0.0464 

0.046350 

3 

12.23 x 10“ 3 

12.2328 x 10 -3 

0.1189 

0.118861 

0.0406 

0.0406266 

4 

12.82 x 10“ 3 

12.81865 x 10 -3 

0.1235 

0.12346 

0.0384 

0.038415 

5 

12.97 x 10“ 3 

12.97 x 10“ 3 

0.1246 

0.124625 

0.0375 

0.03745 

00 

13.02 x 10“ 3 

13.0208 x 10“ 3 

0.1250 

0.1250 

0.0375 

0.0375 


Table. 2: Convergence study for Deflection and Bending moments at the center of simply supported square Kirchhoff plates 

under uniform load 


No of terms 

Wmax 

^xlO- 2 

D 

Mxx max 

pa 2 x 10 -2 

Myy max 

pa 2 x 10 -2 

1 

0.416 

5.34 

5.34 

2 

0.405 

4.69 

4.69 

3 

0.406 

4.86 

4.94 

4 

0.406 

4.81 

4.90 

Exact 

0.406 

4.79 

4.79 
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Table. 3: Simply supported rectangular Kirchhoff plates under point load at the center 


b/a 

Wmax 

Wmai 


Timoshenko and Woinowsky-Krieger 

Present study 

1.0 

0.01160 

0.01160 

1.2 

0.01353 

0.01353 

1.4 

0.01464 

0.01464 

1.6 

0.01570 

0.01570 

1.8 

0.01620 

0.01620 

2 

0.01651 

0.01651 


V. CONCLUSIONS 

The double finite Fourier sine transform method has been 
used to derive analytic flexural solutions of simply 
supported rectangular Kirchhoff plates under general 
distributed load, point load at r|, uniform patch load 
and uniform load on the entire plate. The analysis is 
performed without any assumption of the displacement 
trial (shape) function which illustrates the advantage of 
the method. The method is an efficient and accurate 
analytical tool for Kirchhoff plate bending analysis and 
can be extended to other boundary value problems of 
plates such as buckling and vibration. 
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